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ABSTRACT 


The  deformation  mechanisms  of  aluminum  matrix  composites  containing  Si: 
whiskers  cr  particulate  were  investigated  both  experimentally  and  analytically 
to  determine  factors  critical  to  their  mechanical  behavior.  The  effects  of 
the  discontinuous  SiC  on  age-hardening  behavior;  the  influence  of  SiC 
morphology;  interactions  between  SiC  dispersoids  and  matrix  precipitates;  and 
interactions  between  dispersoids,  matrix  precipitates,  and  dislocations  were 
considered.  Two  matrix  alloys  were  investigated:  a  solution-hardened  5456 
alloy  containing  no  matrix  precipitates;  and  a  2124  alloy  that  could  be  aged 
and  cold  worked  to  contain  different  matrix  precipitate  mi crostructur es  and 
dislocation  contents.  Materials  containing  0$,  8?,  and  20?  by  volume  of 
whisker  or  particulate  were  prepared  for  each  alloy,  using  the  best  powder 
materials  available. 

The  influence  of  SiC  on  the  age-hardening  behavior  of  2124  was  evaluated 
using  differential  scanning  calorimetry.  The  precipitation  sequence  in  the 
composites  was  found  to  be  very  similar  to  that  in  SiC-free  material  prepared 
by  either  powder  or  ingot  metallurgy.  Both  whiskers  and  particulate  increased 
the  quench  sensitivity  of  2124  by  causing  the  precipitation  of  GPB  zones  and 
S’  during  the  quench  from  solution  heat  treatment.  This  effect  was 
independent  of  quench  rate. 

The  tensile  stress-strain  behavior  of  these  composites  was  experimentally 
determined.  Comparisons  were  made  between  the  5456  and  2124  containing  SiC 
whisker  or  particulate,  and  between  2124  materials  having  systematic  changes 
in  matrix  precipitate  and  dislocation  content.  The  elastic  modulus  and  work¬ 
hardening  rate  increased  with  SiC  content.  Whiskers  increased  the  modulus 
more  than  particulate.  The  proportional  limit  of  the  composites  was  generally 
less  than  that  of  the  unreinforced  alioys.  This  is  thought  to  be  due  to 
localized  stress  concentrations  caused  by  the  nondeforming  SiC.  The  state  of 
matrix  precipitation  was  found  to  have  a  dramatic  effect  on  the  proportional 
limit  and  work-hardening  rate  of  the  composite.  The  proportional  limit  for 
the  20?  SiC/2124  increased  from  9  ksi  in  the  annealed  condition  to  95  ksi  for 
the  T8  temper.  These  effects  are  related  to  mechanical  interactions  between 
SiC  dispersoids,  precipitates,  and  dislocations,  and  to  the  effects  of 
dislocations  and  SiC  on  the  precipitate  distribution. 

Two  finite  element  models  (FEM)  utilizing  elastic-plastic  analysis  were 
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constructed  to  describe  the  stress-strain  behavior  of  the  composites.  In 
both,  the  SiC  fibers  are  represented  as  cylinders  with  length  1  and  diameter  d 
in  a  three-dimensional  array.  The  ends  are  aligned  in  one  mcdei,  whereas  they 
are  staggered  in  the  other.  Good  agreement  between  the  calculated  and 
experimental  stress-strain  behavior  of  the  composite  was  obtained  using  the 
experimentally  measured  characteristics  of  the  unreinforced  matrix  alloy  as  a 
starting  point.  Excellent  correlation  was  also  obtained  for  the  effects  of 
1/d  ratio  for  both  particulate  and  whisker  material.  The  FEM  accurately 
predicts  the  observed  increase  in  elastic  modulus  and  decrease  in  proportional 
limit  with  SiC  additions,  as  well  as  the  increase  in  modulus  with  increasing 
1/d  ratio.  The  first  two  of  these  trends  are  incorrectly  predicted  by  shear 
lag  theory.  In  addition,  the  location  of  initial  yielding  was  determined  for 
both  longitudinal  and  transverse  loading  and  was  found  to  be  a  sensitive 
function  of  1/d.  These  models  offer  great  potential  for  predicting  and 
improving  the  mechanical  behavior  of  dis continuously  reinforced  composites. 
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1. OVERALL  SUMMARY  OF  RESULTS 

Aluminum  alloys  reinforoed  with  whisker  or  partiouiate  SiC  have  generated 
great  interest  because  of  their  increased  stiffness  and  potential  for 
increased  strength,  when  compared  with  current  aluminum  alloys.  In  contrast 
to  continuous  filament  composites,  they  also  offer  isotropic  properties, 
easier  fabricability  and  formability,  and  potentially  lower  cost.  However, 
basic  questions  exist  about  the  fundamental  behavior  of  these  materials,  and 
their  properties  are  far  from  optimal.  Since  these  materials  are 
discontinuously  reinforced  and  generally  much  weaker  than  continuous  filament 
material,  the  mechanical  properties  of  the  matrix  alloy  play  a  significant 
role  in  the  overall  properties  of  the  composite.  Currently  obtainable 
strength  levels  are  not  impressive,  so  increased  strength  from  age  hardening 
is  desirable.  In  addition,  discontinuously  reinforced  materials  tend  to  fail 
in  a  relatively  brittle  manner,  and  their  ductility  and  fracture  toughness  are 
low.  It  is  clear  that  the  strength,  ductility,  and  toughness  of  these 
materials  are  controlled  by  interactions  between  the  SiC  particles,  the  age¬ 
hardening  precipitates,  and  deformation  within  the  aluminum  alloy  matrix. 

The  work  reported  herein  has  addressed  these  critical  interactions  by 
evaluating  the  effects  of  SiC  particles  on  precipitation;  and  by  determining 
the  influence  of  SiC  particles,  precipitate  structure,  and  dislocation  content 
on  mechanical  behavior.  Both  experimental  and  analytic  approaches  have  been 
utilized,  and  factors  of  primary  significance  to  elastic  and  plastic  behavior 
have  been  identified.  In  addition,  the  finite  element  models  that  have  been 
produced  are  consistent  with  experimental  observations  and  provide  insights 
that  should  lead  to  improvements  in  the  mechanical  behavior  of  these 
materials . 

1.1  GENERAL  APPROACH 

The  objective  of  this  work  was  to  obtain  a  basic  understanding  of  the 
deformation  mechanisms  that  operate  in  discontinuously  reinforced, 
pr  ->cipi  tat  ion-hardened  aluminum  alloys,  so  that  the  mechanical  behavior  of 
these  materials  can  be  improved.  This  class  of  materials  contains  SiC 
whiskers  (SiCw)  or  SiC  particulate  (SiCp)  and  age-hardening  precipitates.  In 
this  work,  a  systematic  investigation  was  conducted  using  both  experimental 
and  analytic  approaches  to  determine  the  effect  of  the  dispersoid  and  matrix 
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micr  ostructure  on  mechanical  behavior.  Two  systems  were  considered:  a 
solution-hardened  5-456  aluminum  aiioy  containing  no  matrix  precipitates;  and  a 
2124  aluminum  alloy  that  could  be  aged  and  cold  worked  to  contain  different 
matrix  precipitate  and  dislocation  structures.  Materials  containing  0?,  35, 
and  205  by  volume  of  whisker  or  particulate  SiC  were  commercially  prepared  for 
each  alloy  using  the  best  powder  materials  available.  Initially,  the  presence 
of  SiCw  and  SiCp  on  the  precipitation  behavior  of  PM  2124  was  established 
utilizing  the  technique  of  differential  scanning  calorimetry  (DSC)  to 
determine  the  effects  of  quenching  rate  and  aging  practice.  Once  the 
precipitation  characteristics  associated  with  aging  heat  treatments  were 
established,  the  stress-strain  characteri sti  cs  of  PM  2124  material  containing 
different  precipitate-dislocation  structures  and  SiC  contents  were  determined 
and  compared  with  similar  SiC  content  PM  5456  material.  These  experimentally 
determined  relationships  were  also  compared  with  the  results  of  independently 
developed  finite  element  models  of  expected  stress-strain  behavior  and 
initiation  of  localized  plasticity.  Based  on  these  comparisons,  working 
models  have  evolved  that  will  be  used  to  improve  the  properties  of  these 
complex  but  potentially  significant  materials. 

1.2  EFFECT  OF  SiC  ON  THE  PRECIPITATION  BEHAVIOR  OF  PM  2124  ALUMINUM  ALLOY 

Since  the  presence  of  SiC  whiskers  or  partioulate  could  affect 
precipitation  during  the  heat  treatment  of  2124  aluminum  alloy,  it  was 
necessary  to  characterize  the  precipitate  microstructure  resulting  from 
different  quenching  and  aging  treatments.  For  this  task,  DSC  was  used  to 
evaluate  these  effects  quantitatively.  Five  PM  2124  alloys,  i.e.,  containing 
05,  85  or  205  SiCw  or  SiCp,  as  well  as  an  ingot  metallurgy  (IM)  2124  alloy 
baseline  containing  05  SiC  were  examined  to  enable  consider ation  of 
reinforcement  type  and  concentration.  The  heat  treatments  all  began  with  a 
solution  treatment  of  one  hour  at  520°C.  This  temperature  was  chosen  after 
DSC  analysis  of  the  as-received  composites  showed  that  520°C  was  above  the 
last  precipitate  dissolution  peak  and  below  the  onset  of  melting  for  these 
materials.  Three  cooling  rates  were  used  in  order  to  establish  the  severity 
of  quench  sensitivity  in  2124.  The  quenching  media  were  ice  brine  (IB), 
boiling  water  (BW),  and  stiLl  air  (AC).  The  cooling  rates  associated  with 
these  media  were  1400  K/s,  18.5  K/s,  and  2.5  K/s,  respectively.  Four  aging 
conditions  were  chosen  to  correspond  to  the  as-quenched  (W)  temper,  the  room 
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temperature  aged  (T4)  temper,  the  elevated  temperature  aged  (To)  temper,  and 
the  annealed  (0)  temper. 

The  DSC  behavior  of  these  materials  can  generally  be  described  as  an 
initial  exothermic  (formation)  event  in  the  50°C-125°C  temperature  interval. 
This  is  due  to  the  precipitation  of  Gunier  Preston  GPB  zones  from  the 
supersaturated  solid  solution  retained  by  a  rapid  quench  and  was  not  observed 
in  many  of  the  conditions  examined.  The  second  event  is  an  endothermic 
(dissolution)  peak  in  the  temperature  range  150°C  to  240°C.  This  is  due  to 
the  dissolution  of  GPB  zones.  Immediately  following  is  a  formation  peak 
between  240°C  and  300°C  which  is  due  to  the  precipitation  of  S'. 

Subsequently,  over  the  temperature  range  from  300°C  to  500°C  there  is  a  large 
dissolution  peak.  This  is  assumed  to  be  due  to  the  dissolution  of  S'. 

Based  on  the  reactions  and  peak  areas  observed,  the  results  for  the 
different  conditions  examined  indicate  that  the  precipitation  sequence  in  the 
composites  is  very  similar  to  that  of  SiC-free  material;  additionally,  the 
sequences  were  similar  in  SiC-free  2124  prepared  either  by  PM  or  IM.  The 
quench  sensitivity  of  2124  was  affected  by  the  presence  of  SiC.  Both  SiC 
particulate  and  whiskers  increased  the  quench  sensitivity  of  2124  somewhat, 
mainly  by  causing  the  precipitation  of  GPB  zones  and  S'  during  the  quench. 

The  most  pronounced  effect  of  SiC  was  found  to  be  a  decrease  in  the  volume 
fraction  of  GPB  zones  formed.  This  effect  was  independent  of  quenching  rate, 
i.e.,  it  occurred  in  all  conditions.  In  addition,  the  PM  process  was  found  to 
cause  a  small  but  significant  increase  in  the  rate  of  precipitation  of  GPB 
zones  and  S ' . 

1.3  EFFECT  OF  SiC  ON  STRESS-STRAIN  BEHAVIOR 

The  tensile  stress-strain  behavior  of  a  series  of  discontinuously 
reinforced  aluminum  alloy  composites  containing  0%,  8%,  and  20%  by  volume  SiC 
was  examined.  Comparisons  were  made  between  the  behavior  of  age-hardenable 
(2124)  and  solution-hardened  (5456)  matrix  alloys,  whisker  and  particulate 
forms  of  the  SiC  disporsoid  reinforcement,  and  various  states  of  matrix 
precipitation.  The  elastic  modulus  and  work-hardening  rate  of  both  matrix 
alloys  increased  systematically  with  the  addition  of  SiC.  The  effects  of 
whiskers  on  the  modulus  were  greater  than  those  of  particulate.  In  both 
alioys,  the  proportionaL  limit  of  the  8%  SiC  composites  was  less  than  that  of 
the  unreinforced  PM  matrix.  Increasing  the  SiC  to  20%  raised  the  proportional 
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limit,  but  in  2124,  the  proportional  limit  of  the  20%  composites  was  stiii 
below  that  of  the  FM  material.  This  reduction  in  prcpcrticral  limit  beha vie-' 
was  characteristic  of  ail  of  the  heat  treatment  conditions  examined,  vicu  tr.e 
exception  of  the  T8  temper,  which  showed  a  monotonic  increase  in  the 
proportional  limit  with  increased  SiC  content.  The  reduction  in  proportional 
limit  may  be  caused  by  the  presence  of  mobile  dislocations  generated  by  the 
thermal  expansion  mismatch  between  SiC  and  aluminum,  or  by  localized  stress 
concentrations  associated  with  the  SiC  dispersoids.  The  FEM  results  described 
below  indicate  that  the  latter  effect  predominates. 

The  state  of  matrix  precipitation  was  found  to  have  a  pronounced  effect 
on  the  mechanical  properties  of  the  composite.  The  proportional  limit  of  the 
20%  SiC/2124  composite  varied  from  9  ksi  in  the  annealed  condition  to  95  ksi 
in  the  T8.  For  8?  SiC/2124  material,  the  range  was  from  9  to  70  ksi.  The 
other  deformation  characteristics  of  these  materials  (work-hardening  rate, 
yield  stress,  etc)  also  varied,  but  not  so  widely.  The  overall  mechanical 
behavior  of  these  composites  is  understandable  in  terms  of  the  interactions 
between  the  high  density  of  mobile  dislocations,  the  plastic  nonhomegenei ty 
introduced  by  the  nondef ormable  SiC,  the  effects  of  various  precipitate 
distributions  on  dislocation  movement  in  the  matrix,  and  the  effects  of 
dislocations  on  the  precipitate  distribution. 

1 . 4  FINITE  ELEMENT  MODELING 

Our  final  year’s  work  was  concerned  with  creating  and  qualifying  FEMs  of 
short  fiber  reinforced  composites  and  using  these  models  to  predict  mechanical 
behavior.  The  initial  qualification  and  testing  of  the  model  showed  that  the 
stress-strain  characteristics  of  the  composites  predicted  by  the  model  using 
the  measured  properties  of  the  matrix  and  the  SiC  as  input  closely  matched  the 
experimental  observations.  More  detailed  comparisons  of  the  effects  of  SiC 
1/d  ratio  and  volume  fraction  on  the  modulus  and  proportional  limit  of  the 
composite  showed  that  the  model  was  extremely  accurate  in  predicting  both  the 
trends  of  the  data  and  the  actual  numerical  values.  These  results  were  also 
compared  to  the  predictions  of  the  best  available  shear  lag  theories.  It  was 
found  that  the  shear  lag  theories  were  not  even  able  to  predict  accurately  the 
trends  of  the  data,  much  less  the  actual  values.  In  particular,  the  FEM 
accurately  predicted  an  initial  increase  in  the  modulus  of  the  composite  when 
SiC  whiskers  or  particulate  are  added.  The  shear  lag  theories  do  not.  In 
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addition,  the  FEM  accurately  predicted  an  initial  decrease  in  the  proportional 
limit  of  the  composite  when  SiC  whiskers  or  particulate  are  added.  Again,  the 
shear  lag  theories  do  not.  The  behavior  of  the  proportional  limit  was 
predicted  to  be  a  moderately  complex  function  of  the  1/d  ratio  and  the  volume 
fraction  SiC,  and  is  the  result  of  local  constraints  on  deformation  in  the 
matrix  in  the  vicinity  of  the  nondeforming  SiC.  These  predictions  are  unique 
to  the  FEM  model  and  provide  important  new  insights  into  the  mechanical 
behavior  of  these  systems. 

The  model  also  provides  additional  information  on  the  micromechanics  of 
composite  deformation.  Specifically,  the  model  shows  where  in  the  composite 
the  yield  stress  of  the  matrix  is  first  exceeded  during  a  s tress-strain 
test.  This  locus  of  yielding  was  found  to  be  a  sensitive  function  of  the  1/d 
ratio  of  the  reinforcement  when  the  applied  stress  was  in  a  longitudinal 
direction.  The  locus  moved  from  the  midpoint  between  the  particles  at  l/d=1 
to  the  fiber  ends  themselves  when  l/d-2!.  For  transverse  loading,  the  shift  in 
yield  locus  was  not  as  significant.  These  results  are  described  in  greater 
detail  in  the  following  section. 
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2.  SUMMARY  FOR  FINAL  YEAR 


Work  during  the  last  year  was  concerned  with  creating  and  qualifying  FEMs 
of  short  fiber  reinforced  composites,  and  using  these  models  to  predict  the 
behavior  of  various  hypothetical  systems.  Our  objective  in  this  work  was  to 
utilize  the  models  to  help  in  understanding  our  previous  year's  experimental 
data,  and  then,  if  the  models  proved  to  be  reliable,  to  use  them  for 
exploration  of  the  effects  of  various  system  characteristics  (i.e.,  yield 
strength  and  work-hardening  rate  of  the  matrix  alloy,  etc)  on  the  mechanical 
behavior  of  these  composites.  We  were  also  interested  in  using  the  models  to 
predict  how  deformation  would  take  place  on  a  local  scale  iu  the  composite. 
This  information  would  then  be  useful  in  experimental  investigations  of 
deformed  microstructures.  As  shown  below,  the  models  we  developed  correlate 
extremely  well  with  the  experimental  data  and  have  been  used  to  make 
insightful  predictions. 

2.1  FINITE  ELEMENT  MODELS 

Two  FEMs  were  constructed  similar  in  nature  to  the  models  of  previous 
researchers  (Ref  1-2).  In  both  of  the  models,  the  SiC  fibers  are  represented 
as  cylinders  with  a  length  1  and  a  diameter  d.  These  cylinders  are  embedded 
in  a  matrix  volume  element  that  had  a  square  cross  section.  The  transverse 
fiber  spacing  was  D  and  the  longitudinal  fiber  spacing  was  L.  A  sketch  of  the 
first  model  is  shown  in  Fig.  1.  In  this  model,  the  fibers  are  arranged  on  a 
three-dimensional  rectangular  grid,  and  the  ends  are  aligned.  Thus,  the 
volume  element  sketched  in  Fig.  1  represents  one-eighth  of  the  total  volume 
element.  This  is  all  that  is  required  for  our  calculations,  because  the  rest 
of  the  information  is  obtained  by  symmetry.  With  this  volume  element,  both 
longitudinal  and  transverse  properties  can  be  calculated.  An  earlier  model 
was  similar  to  this  one  but  had  a  one-sixteenth  unit.  This  was  adequate  for 
longitudinal  properties,  but  was  not  suitable  for  calculation  of  transverse 
properties;  thus,  it  was  extended  to  the  one-eighth  model  of  Fig.  1. 

The  volume  element  of  the  second  model  is  shown  in  Fig.  2.  This  model  is 
designed  to  arrange  the  fibers  in  a  staggered  manner  so  that  their  ends  could 
overlap  as  shown  in  Fig.  2,  and  was  developed  specifically  to  evaluate  the 
effect  of  a  less  ordered  fiber  arrangement  on  the  resultant  composite 
properties.  Results  from  these  calculations  are  presented  in  Section  2.6. 
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Fig.  2  Volume  Element  of  the  Staggered  Fiber  FEM  Model;  the  Fiber  I/d  ii  Equal  to  the  Matrix  L/D,  but  the  Fibers 
Are  Arranged  Such  That  Their  Ends  Are  Staggered. 
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Ali  of  cur  models  were  constructed  on  the  PATRAN  graphics  system  (Ref 
3).  An  inhouse  translator  was  then  used  to  convert  the  PATRAN  output  into  an 
input  data  deck,  including  the  multipoint  constraint  equations  along  the 
planes  of  symmetry,  for  use  with  our  three-dimensional  finite-element  program 
"HES"  (Ref  JO.  This  program  was  then  used  to  perform  elastic-plastic  analyses 
on  the  constructed  models.  The  output  of  the  computer  runs  included  a 
uniaxial  tensile  stress-strain  curve  and  maps  of  the  strains  and  stresses  in 
each  element.  The  input  data  included  the  stress-strain  characteristics  of 
the  matrix  alloy,  which  were  obtained  from  our  experimental  data  on  the 
identically  processed  SiC-free  materials,  and  the  elastic  modulus  of  the  SiC 
whiskers/particulates,  which  was  assumed  to  be  485  GPa.  The  results  of  these 
calculations  are  described  below. 

2.2  EFFECTS  OF  Pc  AND  ct  ON  CALCULATED  STRESS-STRAIN  CURVES 

After  the  model  was  generated  and  validated,  the  first  task  undertaken 
was  an  examination  of  the  sensitivity  of  the  predicted  composite  stress-strain 
curve  to  the  deformation  characteristics  of  the  matrix  alloy.  The 
experimentally  measured  stress-strain  curve  of  the  2124-T4  PM  matrix  material 
was  approximated  in  the  model  as  a  bilinear  curve  with  initial  slope  (modulus) 
E,  proportional  limit,  ?Q,  and  final  slope,  ET,  the  work-hardening  rate.  For 
computational  convenience,  the  final  slope  was  input  as  a  ratio  a,  where  a  is 
defined  as  E^  divided  by  E,  i.e.,  a»Ej/E.  The  measured  values  of  these 
parameters  were:  E*73  GPa,  Pc»385  MPa;  and  a-0. 1 .  Some  of  the  calculated 
stress-strain  curves  for  l/d-4  are  plotted  in  Fig.  3,  along  with  the 
experimentally  observed  stress-strain  curves  for  PM  2124/-T4,  a  20$  SiCw/2124- 
T4  composite,  and  the  assumed  SiC  behavior.  The  upper  three  calculated 
stress-strain  curves  show  the  sensitivity  of  the  model  to  a,  the  work- 
hardening  rate.  These  three  curves  have  matrix  properties  of  Pc»385  MPa  and 
ot»0  (perfectly  plastic),  a-0.1  or  a»0.2.  It  can  be  seen  that  increasing  the 
work-hardening  rate  of  the  matrix  also  increases  the  work-hardening  rate  of 
the  composite  while  leaving  the  proportional  limit  unchanged.  The  next  two 
curves  show  the  effect  of  keeping  o  constant  at  0.1  while  varying  P.,.  In  this 
case,  a  decreased  matrix  proportional  limit  results  in  a  lower  composite  work- 
hardening  rate  and  a  lower  composite  proportional  limit.  In  ali  cases,  the 
calculated  composite  elastic  modulus  agreed  well  with  the  observed  value  and 
was  unaffected  by  these  parameters.  In  order  to  match  exactly  the  observed 


STRAIN  1%) 


R87-44S3-003(T) 

Fig.  3  Comparison  Between  the  Experimental  Stress-Strain  Curves  for  2124-T4  and  20% 
SiC/21 24-T4,  and  Stress-Strain  Curves  Calculated  by  FEM  Showing  the  Effects  of 
Varying  the  Proportional  Limit,  Pc,  and  Work  Hardening  Rate,  a,  of  the  Matrix. 
{I/d  *  4.)  The  calculated  curve  with  Pc  »  289  and  a  *  0.05  is  an  exact  match  to 
the  experimentally  obtained  20%  SiCw/2124-T4  curve. 


20$  SiCw/2124-T4  composite  stress-strain  curve,  the  matrix  characteristics  had 
to  be  set  at  P„»289  MPa  and  a=0.05.  This  curve  is  also  shown  in  Fig  3  and 
provides  the  best  fit  to  the  experimental  observations  f or  a  unidirectional 
model  with  l/d=4.  These  matrix  characteristics  are  lower  than  those  observed 
in  the  SiC-free,  control  material.  A  possible  explanation  is  that  the 
presence  of  SiC  whiskers  in  the  composite  alters  the  stress-strain 
characteristics  of  the  matrix  alloy  (i.e.,  dislocation  density,  precipitate 
distribution,  residual  stress,  eto).  It  is  well  established  from  other 
experiments  that  alteration  to  the  matrix  microstructure  takes  place;  thus,  it 
appears  that  the  matrix  in  the  composite  has  a  lower  flow  stress  and  a  lower 
work-hardening  rate  than  found  in  the  SiC-free  PM  material. 

2.3  EFFECTS  OF  1/d  ON  CALCULATED  STRESS-STRAIN  CURVES 

The  effeccs  of  varying  1/d  between  1  and  10  are  shown  in  Fig.  4  for  a  20$ 
SiC^/5456  composite.  The  experimental  data  for  PM  5456,  20$  SiCw/5456,  and 
20$  SiCp/5456  are  also  shown  in  Fig.  4  along  with  the  assumed  SiC  curve.  It 
can  be  seen  that  increasing  the  1/d  ratio  increases  the  modulus,  proportional 
limit,  and  work-hardening  rate  of  the  calculated  stress-strain  curves.  The 
calculated  stress-strain  curve  for  1/d -1  is  very  similar  to  the  experimental 
curve  from  the  particulate-loaded  material.  The  major  region  of  disagreement 
between  the  two  curves  lies  in  the  initial  yield  region,  where  the 
experimental  data  show  a  more  rounded  behavior.  This  rounded  yield  behavior 
may  be  due  to  the  introduction  of  a  large  number  of  mobile  dislocations  by  the 
differential  thermal  contraction  of  the  SiC  and  aluminum.  Thus,  the  matrix 
material  in  the  composite  is  likely  to  have  somewhat  different  deformation 
characteristics  from  the  matrix  material  without  SiC  that  is  used  as  our 
model,  leading  to  the  discrepancy  in  this  region.  Overall,  however,  the 
agreement  between  the  calculation  and  the  experiment  is  excellent. 

The  results  calculated  for  1/d  greater  than  1  are  also  shown  in  Fig.  4, 
along  with  the  experimental  data  from  the  whisker-loaded  composite.  For  the 
20$  SiCw/5456,  the  best  agreement  appears  to  be  for  an  assumed  l/d-2  or  3-  In 
order  to  assess  the  validity  of  this  assumption,  the  whisker  1/d  was 
measured.  Whiskers  were  extracted  by  an  acid  digestion  of  the  matrix  and 
examined  in  a  scanning  electron  microscope.  The  observed  distribution  of  1  /d 
values  was  found  to  have  a  mean  value  of  4.1  and  a  most  probable  value  of 
2.5.  Thus,  the  agreement  between  the  FEM  predictions  for  1/d  and  the  data  is 
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Fig.  4  Comparison  Between  the  Experimental  Stress-Strain  Curves  for  5456-W,  20%  SiCp/5456-W  and 
SiC^/5456-W.  and  FEM  Calculations  Showing  the  Effects  of  I/d  on  the  Results. 


excellent . 


The  FEM  calculations  provide  a  prediction  of  the  modulus  and  proportional 
limit  of  the  composite.  It  is  interesting  to  compare  these  predictions  to  the 
measured  values  and  to  other  predictions.  Calculated  values  of  the  modulus 
and  proportional  limit  are  listed  in  Table  1  for  both  8?  and  20?  SiC,  along 
with  the  experimentally  observed  values.  In  addition,  Table  1  also  lists  the 
modulus  and  proportional  limit  calculated  using  shear  lag  theory.  For  a 
description  of  these  theories,  see  Ref  5.  The  FEM  calculations  predict  an 
increase  in  the  modulus  of  the  composite  with  additions  of  SiC  and  with 
increasing  1/d  ratios.  The  predicted  values  are  in  excellent  agreement  with 
the  experimentally  measured  ones.  In  contrast,  the  shear  lag  theory 
calculations,  using  either  of  two  formulations,  predict  that  the  modulus  will 
decrease  when  particulate  or  short  fibers  are  added  to  the  matrix  and  will 
increase  only  when  the  1/d  ratio  of  the  SiC  exceeds  approximately  3.  This 
prediction  is  contrary  to  experimental  observations  and  is  due  to  the  shear 
lag  assumption  that  the  fiber  ends  are  debonded. 


The  FEM  calculations  predict  that  the  proportional  limit  will  decrease 
when  SiC  is  added,  primarily  due  to  localized  stress  concentrations.  This  is 
in  agreement  with  the  experimental  observations.  In  contrast,  the  modified 
shear  lag  theory  predicts  that  the  proportional  limit  will  increase  when  SiC 
is  added.  The  shear  lag  theory  predictions  are  significantly  higher  than  the 
experimental  observations  and  the  FEM  results.  Thus,  these  comparisons  show 
that,  for  both  the  modulus  and  proportional  limit,  the  FEM  calculations 
accurately  predict  the  observed  behavior,  while  shear  lag  calculations  are 
quantitatively  and  qualitatively  incorrect. 

2.4  RESULTS  FOR  L/D  *  1/d 


An  assumption  made  for  all  of  the  preceding  calculations  was  that  L/D  - 
1/d  for  the  model  shown  in  Fig.  1.  This  fixes  the  geometry  of  the  element  and 
controls  the  f  iber-to-f iber  end  spacings.  In  order  to  examine  the 
implications  of  this  assumption  on  the  results,  several  calculations  were 
performed  for  different  fiber-end  spacings.  The  results  of  some  of  these  are 
shown  in  Fig.  5,  where  the  results  from  a  calculation  in  which  1/d  »  4  and  L/D 
»  4  or  2.3  are  plotted  along  with  the  experimental  data  for  the  20?  SiC/5456 
composites.  (When  L/D  -  2.3,  the  separation  between  the  fiber  ends  was  0.6d; 
when  L/D  *  4,  the  separation  was  2.3d.)  As  can  be  seen  in  Fig.  5,  differences 
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Fig.  5  The  Effect  of  Changing  the  Assumed  L/D  while  Keeping  I/d  ■  4;  when  L/D  ■  4,  the  Distance  Between 
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between  the  two  predictions  are  evident,  but  they  are  not  particularly 
significant  in  this  case. 

2.5  STRESS-STRAIN  CURVES  FOR  TRANSVERSE  LOADING 

Although  the  whiskers  in  these  composites  are  initially  random,  the 
extrusion  process  causes  significant  alignment  in  the  extrusion  direction. 
Thus,  the  composites  exhibit  a  preferred  whisker  alignment  along  the 
longitudinal  direction,  but  this  alignment  is  not  perfect.  All  of  the 
calculations  described  earlier  were  made  for  a  perfectly  aligned  composite 
with  the  stress  in  the  longitudinal  direction;  thus,  it  is  of  interest  to 
examine  the  tensile  behavior  in  the  transverse  direction,  both  for  comparison 
to  longitudinal  behavior  and  for  an  indication  of  hew  imperfectly  aligned 
whiskers  would  affect  the  results.  The  calculated  transverse  str ess-strain 
curve  is  shown  in  Fig.  6  for  20%  SiC/5456  with  1/d  ratios  between  1  and  10. 

The  results  show  that  there  is  no  effect  of  1/d  on  the  transverse  stress- 
strain  behavior.  All  of  the  curves  are  virtually  identical.  Comparison  of 
these  results  with  the  calculated  longitudinal  results  in  Fig.  4  shows  that 
the  calculated  transverse  stress-strain  curve  lies  below  the  longitudinal 
results  and  that  the  longitudinal  result  is  much  closer  to  the  experimental 
data.  This  is  consistent  with  the  observation  that  there  is  a  significant 
amount  of  fiber  alignment,  and  also  indicates  that  the  effect  of  misalignment 
will  be  to  lower  the  stress-strain  curve. 

2.6  RESULTS  FROM  THE  STAGGERED-FIBER  MODEL 

The  calculations  described  earlier  all  were  performed  with  the  first 
model  where  the  fiber  ends  are  aligned.  Calculations  were  also  made  using  the 
staggered-f iber  model  shown  in  Fig.  2.  Results  from  these  calculations  for  a 
20%  SiC/5456  composite  with  l/d-1,  and  10  are  shown  in  Fig.  7.  A 
comparison  of  these  stress-strain  curves  to  those  obtained  with  the  aligned 
fiber  model  (Fig  4)  shows  that  the  aligned  fiber  model  gives  a  better  overall 
fit  to  the  particulate  (l/d-1)  data,  but  that  the  staggered-f iber  model  fits 
better  in  the  initial  yield  region.  This  indicates  that  the  local  constraints 
imposed  by  fiber  overlap  may  play  a  significant  role  in  yielding  behavior. 
There  are  two  differences  between  the  aligned  and  staggered  models.  First, 
the  staggered  arrangement  has  the  fibers  more  closely  packed,  the  separation 
between  two  adjacent  fibers  in  the  transverse  direction  being  /  2  closer,  and 
secondly,  the  fiber  ends  overlap  for  higher  1/d  ratios.  For  high  1/d  ratios 
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Fig*  6  Calculated  Stress-Strain  Curves  for  Transverse  Loading  of  a  20%  SiC/5456  Composite  with  Various  I/d 
Ratios;  the  Experimental  Data  Are  for  Longitudinal  Loading. 
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Fig.  7  Calculated  Stress-Strain  Curves  for  Longitudinal  Loading  of  the  Staggered  Fiber  Model  Compared  to 
Experimental  Data  From  5456. 


(i/d  greater  than  approximately  *0),  the  overlapping  effect  causes  a  'fiber  tip 
to  "see"  the  center  region  of  an  adjacent  fiber,  and,  thus,  there  is  a 
tendency  to  reduce  the  peak  stress  around  the  fiber  tip.  For  iow  I'd  ratios, 
the  fiber  tips  tend  to  "see"  one  another,  thus  raising  the  stress  level  around 

the  fiber  tip.  As  can  be  seen  in  Fig.  4  and  7  and  Table  2,  for  low  aspect 

ratios,  the  staggered  model  has  less  work  hardening  and  yields  sooner  than  its 
aligned  model  counterpart,  while  for  l/d-10,  the  staggered  model  has  mere  work 

hardening  and  yields  later.  The  work  of  Ref  2  shows  similar  findings  for  high 

aspect  ratios,  but  low  aspect  ratios  were  not  studied.  The  lower  proportional 

limit  predicted  by  the  staggered-f iber  model  reinforces  our  previous 

observation  that  local  stress  enhancement  due  to  stress  concentrations  at  the 
fiber  tips  is  the  primary  cause  of  the  observed  decrease  in  proportional  iimit 
when  Sic  is  added. 

2.7.  LOCATION  OF  INITIAL  YIELDING 

The  calculations  described  earlier  contain  a  large  amount  of  additional 
information.  One  aspect  of  particular  interest  is  the  location  of  the  volume 
element  in  the  matrix  where  yielding  first  takes  place.  Some  of  these  data 
are  shown  in  Fig.  8  and  9,  where  the  location  of  initial  yielding  is  indicated 
for  the  aligned  fiber  model  in  longitudinal  and  transverse  tension  with  l/d=1 , 
2,  and  4.  For  the  longitudinal  case,  it  can  be  seen  that  the  location  of 

yielding  moves  away  from  the  midpoint  in  the  matrix  between  the  fiber  ends,  to 

the  circular  edges  of  the  fibers  themselves.  This  is  a  significant  resuit  for 
understanding  the  differences  in  behavior  of  particulate-  and  whisker- 
reinforced  composites,  and  indicates  that  the  micromechanical  flow  and 
fracture  mechanisms  are  different  in  the  two  materials.  In  the  transverse 
case,  there  is  also  a  dependence  of  the  location  of  initial  yielding  on  the 
1/d  ratio,  as  shown  in  Fig.  9.  However,  in  the  transversely  loaded  case, 
initial  yielding  always  occurs  between  adjacent  fibers. 


£ 

% 

a 


X?' 


$ 

is 


Table  2  Effect  of  S/d  on  Modulus  (E)  and  Proportional 
Limit  (Pc)  of  20%  SiCd/5456 


I/d 

Aligned  fiber  model 

Staggered  fiber  model  j 

E.  GPa 

Pc  MPa 

E,  GPa 

Pc  MPa 

1 

110 

181 

108 

170 

2 

114 

217 

114 

188 

4 

124 

224 

128 

222 

10 

129 

242 

144 

275 
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